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ABSTRACT 



Aims. We study the soft X-ray emission induced by charge exchange (CX) collisions between solar-wind, highly charged ions and 
neutral atoms of the Martian exosphere. 

Methods. A 3D multi species hybrid simulation model with improved spatial resolution (130 km) is used to describe the interaction 
between the solar wind and the Martian neutrals. We calculated velocity and density distributions of the solar wind plasma in the 
Martian environment with realistic planetary ions description, using spherically symmetric exospheric H and O profiles. Following 
that, a 3D test-particle model was developed to compute the X-ray emission produced by CX collisions between neutrals and solar 
wind minor ions. The model results are compared to XMM-Newton observations of Mars. 

Results. We calculate projected X-ray emission maps for the XMM-Newton observing conditions and demonstrate how the X-ray 
emission reflects the Martian electromagnetic structure in accordance with the observed X-ray images. Our maps confirm that X-ray 
images are a powerful tool for the study of solar wind - planetary interfaces. However, the simulation results reveal several quantitative 
discrepancies compared to the observations. Typical solar wind and neutral coronae conditions corresponding to the 2003 observation 
period of Mars cannot reproduce the high luminosity or the corresponding very extended halo observed with XMM-Newton. Potential 
explanations of these discrepancies are discussed. 
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1. Introduction 



Sever al classes of solar system bodies, comets (L isse et al 



1996I) planets dCravensI l2000b iDennerll |2003: iFuiimoto etal. 
2007 : ISnowden et alj2009TT and even th e interplaneta ry med ium 
dCoxlll998h ISnowden et"aLll2004t iKoutroumpa et al.l l2007h are 



known sources of X-ray radiation. The X-ray emission common 
to all those objects is produced from the radiative de-excitation 
of highly-charged solar wind (SW) ions that have captured an 
electron from neutral particles encountered in these environ- 
ments, and is accordingly named solar wind charge exchange 
(SWCX) X-ray emission. A thor ough review of S W CX emis- 
sion in the solar system is given in Bh ardwajl j2010h and lDennerl 
(I2010h . 

In this paper, we focus on CX-induced X-ray emission 
from nonmagnetized planets, specifically from Mars. Martian 
CX e mission was predicted for the first time by ICravensi 
(l2000l) . and the first estimates of the Martian SWCX emission 
based on comet obs ervations and calculations were given by 
iKrasnopolskvl (l200lh . Following that, the first model of minor 
ion propagation throug h a simplified SW - Mar tian exosphere 
interface was devised bv lHolmstrom et al] J2001 ). Thanks to the 
absence of an internal magnetic field, the SW can interact di- 
rectly with the upper martian atmosphere, and CX collisions be- 
tween multiply charged heavy SW ions and atmospheric neutral 
constituents may occur at high altitudes. The first detection of 



Marti an X-rays was confirmed very soon by Chandra observa- 
tions (lDenner]|l2002l) . In addition to a fully illuminated disk of 
the size of the planet, a fainter halo was detected up to three 
Martian radii. Fluorescence and scattering of solar X rays are 
efficient processes at low altitudes, and they explain emissions 
from the disk, while X-ray inducing CX collisions provide a 
good description of the halo emission. 

iGunell et alj d2004l) estimated the CX X-ray emission from 
Mars in the context of the Chandra observations, incorporating 
a hybrid model for the SW-Mars interaction and a test particle 
simulation of heavy ion trajectories near Mars. They simplified 
the computation by considering a unique cross section for both H 
and O atoms, and replacing cascading transitions to ground state 
by a simplified two-step cascade model. This model assumes 
equal probabilities for all radiative transitions starting from the 
same initial excited state. The simulation shows good agreement 
with the Chandra 2001 observation on the total X-ray intensity. 
Chandra's limited spectral resolution, however, did not allow any 
conclusive spectral analysis. 

Mars was also observed with the Refle ction Grating 
Spect rometer (RGS) aboard XMM-Newton dDennerl et al.l 
120061) . RGS' spectral resolution allowed for the first time a de- 
tailed spectroscopic study of Martian X-rays from the disk and 
halo, resolving fluorescence spectral lines induced in the up- 
per atmosphere from SWCX emission lines produced in colli- 
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sions between SW heavy ions and exospheric neutrals at greater 
distances from the planet. This analysis yielded the strongest 
SWCX emission levels from Mars ever observed, with emission 
being detected up to eight Mars radii, and a total halo luminosity 
of (12.8+1.4) MW (mainly produc e d by S WCX emission lines). 

More recently, Is hikawa et alj d201 lh have observed Mars 
with the Suzaku XIS instrument in 2008, during the minimum 
phase of solar activity. However, their results did not show any 
significant SWCX emission from the planet, and only yielded a 
2<r upper limit of 4.3xl0~ 5 photons crrT 2 s _1 from the oxygen 
lines in the 0.5-0.65 keV energy range. In the same energy range, 
the Chandra and XMM-Newton observations yielded a solid lcr 
detection of 5.8xl0~ 6 photons cirT 2 s _1 and 3.6x1 -5 photons 
cm' 2 s respectively. 

We present an improved model of the SWCX emission in 
Mars' environment, with an unprecedent spatial resolution (130 
km) and detailed spectral information on the SWCX induced 
transitions that allows a direct comparison of the simulation re- 
sults to the XMM-Newton measured photon fluxes and total lu- 
minosities. In section [2] we describe the simulation model and 
initial solar and SW parameters appropriate for the XMM obser- 
vation period derived from near-Earth instruments and extrapo- 
lated to Mars' orbit. In section[3]we present projected emission 
maps, spectral line intensities, and total luminosities resulting 
from the simulations. Finally, in section [4] we discuss how the 
model results compare to the data and previous simulations and 
offer some explanations for the discrepancies we found. 



2. Simulation model 

The simulation process is separated into two main steps. In the 
first we use a three-dimensional (3D) and multi species hybrid 
model to provide a full description of the S W interaction with the 
Martian ne utral environment. This sim ulation model has been 
designed bv lModolo et alJ (l2005ll2006T) : the hybrid kernel of the 
simulation program makes use of the Current Advan ce Method 
and Cy clic Leapfrog scheme (C AM-CL) designed by iMatthewsl 
(1 19941) . We resolve the Maxwell equations and self-consistently 
obtain a quasi-stationary solution to the global electromagnetic 
(EM) field as well as a description of the plasma dynamic in 
the vicinity of the planet. This simulation model successfully 
reproduces the plasma environment surrounding Mars, such as 
the Martian bow shock (hereafter MBS) and the Magnetic pile- 
up boundary (hereafter MPB) obs erved by the Mars G lobal 
Surveyor and Phobos-2 spacecrafts (iTrotignon et al.ll2006h . 

In a second step, once we achieve equilibrium in the hybrid 
simulation, the electric and magnetic field 3D distributions ob- 
tained with the hybrid model are used as the input field for a 
test-particle model. Trajectories of highly charged SW ions are 
followed in the simulation box, and X-ray emissions due to CX 
collisions with the planetary neutrals are computed. A similar 
procedure has been successfully used to investi gate the capture 
of alp ha particles by the Martian atmosphere (IChanteur et alJ 
120091) . 

2. ) . The hybrid model 

The hybrid model has been extensively descri bed in 
iModolo et al.1 J20051 [2006b and iKallio et alJ J20ll . Solar 
and planetary ions are represented by sets of macroparticles, 
describing the full dynamics of each ion species, while electrons 
are described by a massless fluid that ensures the conservation 
of the charge neutrality of the plasma and contributes to the 
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Fig. 1. Relative position of Earth and Mars on 
November 20, 2003, generated by the IMCCE 
(http://www.imcce.fr/en/ephemerides/). A color version of 
this figure is available in the online journal.) 



electronic pressure and the electric currents. The magnetic field 
undergoes a temporal evolution according to Faraday's law. 

The coordinate system is defined such that the X axis points 
away from the Sun, (X = V m ./ HVjJI) is aligned with the local 
SW direction, where V sw is the SW velocity. The Y axis lies in 
the ecliptic, perpendicular to V m and pointing ahead on Mars' 
orbit. The interplanetary magnetic field (IMF) vector lies in the 
XY (ecliptic) plane, following Parker's spiral (</> = -56° with re- 
spect to X axis at Mars). Finally the Z axis is toward the motional 
electric field, defined by E com , = -V sw x Bjmf, and therefore co- 
incides with the north ecliptic pole. 

The spatial resolution of the 3D uniform simulation grid is 
130 km, equal to the inertial length of the SW protons. The 
planet is modeled as a fully absorbing obstacle with a radius 
of 3400 km, and ions penetrating in the obstacle boundary are 
stopped. Open boundaries are used in the SW direction and pe- 
riodic boundaries in the perpendicular directions, except for the 
planetary ions, which are escaping freely from the simulation 
domain. At the first step of the simulation a computation is run 
from time t = to time t = 1000 s, while a nearly stationary 
solution is obtained around time t = 450 s. 

This hybrid simulation is performed for SW and EUV con- 
ditions corresponding to the period of the XMM-Newton obser- 
vation of Mars starting on Nov. 19, 2003 23:50 UT. The differ- 
ent parameters are reviewed below and summarized. The SW 
plasma parameters used for this simulation are provided by the 
Advanced Composition Explorer (ACE) and Wind experiment 
and corrected to the Martian orbit. 

2.1.1. Solar parameters 

The ACE and Wind spacecrafts are located at the Sun-Earth's LI 
point, and Figure[T]presents the relative position for Mars, Earth, 
and Sun at the time of the XMM observation. At first order, Mars 
and the Earth are in the same sector, meaning that near-Earth 
solar and SW monitoring can be used as a proxy to determine 
the SW and EUV conditions near the planet. 
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Fig. 2. Solar parameters measured by near-Earth instruments 
during the three months preceding the XMM-Newton observa- 
tion of Mars in November 2003. From top to bottom the curves 
represent, averaged over 3 solar rotations (81 days): proton ve- 
locity in km s (dot-dashed), the F10.7-cm flux index (dashed), 
proton density in crrT 3 (solid) and +7 /0 relative abundance 
(dotted). In-situ data are from ACE at LI, while the radio flux 
indices were taken from the NOAA Space Weather Prediction 
Center. The solid vertical lines show the start and end of the 
XMM-Newton observation of Mars, and the dashed vertical 
lines show the SW event window applied to the simulations, 
once the Earth-Mars SW delay time is taken into account. 



In Figure|2]we present the F10.7-cm Radio flux index from 
the NOAA Space Weather Prediction CenteiQ and SW in-situ 
data from ACE Science Center Level 2 databas^ for the pe- 
riod of November 2003, when the XMM-Newton observation 
of Mars took place. If we ignore the longitudinal difference in 
Mars' position compared to Earth (ACE) and assume a constant 
SW propagation velocity between the LI point and Mars, then 
the time delay that needs to be applied to the ACE in-situ data 
in order to find the appropriate SW event window at Mars is 
At = Ar/V sw , where Ar = ym - fE is the heliocentric distance 
difference between Mars and the LI point. We find that this 
corresponds approximately to the window between 2003/11/18 
12:00 UT and 2003/11/19 06:00 UT (Figure 0. We derive the 
SW parameters by averaging the ACE data over this window 
and scaling to Mars' heliocentric distance for the particle densi- 
ties (1/r 2 ). 

The main ion species in the SW is H + and simulation values 
have been normalized with respect to H + parameters. The main 
SW plasma parameters adjusted to the Martian orbit in the event 
window defined above are 



Density : 
Bulk speed : 
Thermal speed : 
Electron temperature 
He ++ abundance 
IMF : 



n sw ~ 2 cm 3 
V m ~ 675 km s" 1 
v t h ~ 80 km s _1 
T e = 3. 10 5 K 
4% 

WBimfW = 3 nT, 4> 



= -56° 



The F10.7-cm index can be used as proxy for the EUV so- 
lar conditions that influence the planet's neutral environment. 
During the November 20, 2003 period the F10.7-cm index is 175 
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Fig. 3. Neutral density profiles for the CO2 (magenta), H (red), 
and O (black) coronas used in our simulations. The H2 (blue 
dashed ) and He (green dashed) coronas from iRrasnopolskv 
(|2010[) and He (green dash-dotted) corona from iFox & Had 
(120091) are also plotted for comparison. (A color version of this 
figure is available in the online journal.) 

and the 81 day (three solar rotations) average is 130, suggesting 
that the conditions are close to solar maximum. The neutral en- 
vironment is then chosen adequately. 

2.1.2. Neutral atmosphere and exosphere 

The Martian neutral environment is described by three coronae 
of CO2, H, and O show in Figure|3] Both H2 and He coronas are 
also thought to be important at altitudes from aroun d 400 km to 
1000 km (also presented in Figur e [3] extracted from lFox & Had 
120091: lKrasnopolskvll2002l 120101) . However, the density profiles 
of these neutral populations decrease very quickly with altitude 
and should only have an effect in the X-ray emission in the disk 
region, as we discuss in Section [4] This version of the model 
does not include the H2 and He coronas, but future work will 
include these components. 

All information is based on input from the "Solar Wind 
Model Challenge" international team at the Inter national Space 
Science Institute (ISSI) in Bern (iBrain et al.ll2010l) . Atmospheric 
and exospheric coronas are partly deduced from the Mars ther- 
mosphere global circulation model (MTGCM, iBougher et al.l 
2000, 120061 120081) . unless noted otherwise, assuming a spheri- 
cally symmetric neutral environment. The MTGCM was run for 
Martian longitude position Ls = 270 and F10.7 = 105. Although 
observed EUV conditions are slightly more intense than those 
used for the atmospheric simulations, these simulations are the 
most approriate runs. 

For CO2, the density profile used is 

n COl =5.88 \Q™exp{-j^ + 3.55 10° exp , (D 

which is in very goo d agreement wit h the models of 
iKrasnopolskvl d2002l) and lFox & Hacl (120091) . In the equation, z 
is the altitude in km and n is the density in cm 4 . 

For O, the 'cold' oxygen density profile used represents a fit 
from the MTGCM results: 



1 http://www.swpc.noaa.gov/Data/index.html 

2 http://www.srl.caltech.edu/ACE/ASC/level2/ 



-.cold 



233 1013 e A-Th) + 2M 6XP (-4837)- (2) 
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Below 400 km the computed density profile for thermal O 
atoms agrees well w ith results of the photo chemical model by 
Krasnopolskyl d20 1 Ol) . The 'hot' density profile used is a fit from 
Valeille et aiTd2010l) results: 



Table 1. Parameters of SW minor ions included in simulations 



«- = 1.56 10 4 ^(-^) + 2.92 103 expf-^) 

+ 5 .0110 4 ex P (-^y. (3) 

For hyd rogen, the density profile use d is a combination o f 
results from l Anderson & Hordld 19711) and lKrasnopolskvl (l2002h : 

riH = 10 3 exp 



9.25 10 3 



1 



1 



+3 I0 4 exp 



1.48 10 4 



z + 3393.5 3593.5 
1 1 



z + 3393.5 3593.5 



(4) 



2.2. Test particle simulation 

The second step consists of running a test particle simulation 
for each minor species of the SW with the electric and magnetic 
fields obtained from the hybrid simulation model. We consider 
that the E and B fields are in equilibrium so we no longer intro- 
duce H + and He ++ ions in the simulation box, but only minor 
SW ions (C, N, O). In the test particle simulation for a given mi- 
nor species of the SW, we follow the trajectories of N= 2 x 10 6 
macroparticles having the same initial statistical weight when 
injected into the simulation box. This initial statistical weight 
is computed by considering that the physical flux correspond- 
ing to the N macroparticles injected during one second should 
be equal to the nominal flux of the given minor species in the 
SW. At each time step, meaning at each position of the parti- 
cle, CX collisions are computed in the cell, and the statistical 
weight of the test-particle is reduced to mimic the loss of the 
parent ions caused by collisions. We inject one heavy SW ion 
macroparticle after the other in the box and solve equations of 
motion to calculate its trajectory. The next macroparticle is in- 
jected when the former one has either left the box, is absorbed 
by the planet/obstacle, or its weight is reduced to because of 
successive CX collisions. This technique has been successfully 
applied to investigating the h elium budget in the Martian atmo- 
sphere dChanteur et al.ll2009l) .The reaction the SW ion X +q un- 
dergoes in the simulation is 

X +q + M -> X Mq - l) * + M + (5) 

where M = [H, O] are the major neutral target atoms in the 
Martian environment. In each grid cell we calculate the proba- 
bility of ion X +q capturing an electron from exospheric neutrals, 
register the ion X +<?_1) * production rate, and cumulate the con- 
tributions of ions X +(9_1) * produced by all N X +q ions in the cell. 

As long as it may collide with exospheric neutrals, i.e. in 
practice as long as it remains within the simulation box around 
Mars, the X + ^ ?-1 ' ion produced in reaction|5]can continue to ex- 
change charges until it is neutralized. We investigated the in- 
fluence of these secondary CX collisions for the case of +6 
ions sequentially produced by +8 and +7 ions in Mars' ex- 
osphere. The SWCX emission produced by secondary CX colli- 
sions has been estimated and was found to contribute to the total 
+6 (OVII) emission by less than 1 %, hence can be ignored. 

For all minor SW species the initial distribution function 
in the SW is an isotropic Maxwellian with a bulk velocity of 
675 km/s. Although the SW minor ions do not influence the 



Ion 


[X+'/O] Abundance 


Cross section (10 cm^) 




OlOW o V V 


Real-Time^ 


CX with H CX with O c 




0.318 


0.13±0.06 


4.2 5 


c+ 5 
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0.37+0.03 


4.0 2 


N +7 


0.006 


no data 


5.7 12 


N +6 


0.058 


no data 


3.7 5 


0+8 


0.07 


<0.002 


5.7 6 


o +7 


0.20 


0.04±0.02 


3.4 12 


He ++ 


78±4 


92+3 




He/H 


0.044+0.003 


0.045 





Notes. (< " From ISchwadron & CravensI d2000h {b) From ACE-SWICS 
w From ISchwadrmTlsrc7avmsn2*0*00l) . based on the assumption that 
atomic O has similar weight to HiO. 



EM structure equilibrium around Mars, the SWCX line emis- 
sion produced by the reactions described above is directly pro- 
portional to the parent ion flux. More precisely, the SWCX line 
photon flux is the integral over the line of sight (here limited to 
the simulation box size) of the product of the minor ion density, 
the neutral density (here the exospheric H and O), the relative 
velocity between the two colliding particles (usually averaged to 
the SW ion velocities), the CX collision cross-section, and the 
individua l line emission probability (base d on theoretical calcu- 
lations by [Kharchenko & Dalgarno1l2000h . 

In Table Q] we summarize the parameters for the minor S W 
ions included in the simulations. The first column lists the ion 
charge state. Columns 2 and 3 list the charge state relative abun- 
dances with respect to ox ygen, and alpha (He ++ ) to pro ton ratio 
for the typical slow SW ( ISchwadron & CravensI l200"ol) and the 
real-time values extracted from the ACE database for the win- 
dow of the XMM observation of Mars, as described in section 
12.1.11 The ion density N, is then defined as the product: 



Ni = 



X +q 



O 



0_ 

We 



He 



H 



(6) 



In columns 4 and 5 of Table [T] we also summarize the cross 
sections used in the simulations for CX reactions with H and 
O atoms, respectively. The cross section va lues for H are as- 
sembl ed from references listed in Table 1 of IKoutroumpa et ail 
(120061) . For the cro ss sections in CX with O ato ms, we adopt 
the values listed in ISchwadron & CravensI (120001) for reactions 
with H2O molecules. All cross sections refer to typical slow 
SW velocities of 450 km s _1 . The CX cross sections are known 
to vary smoo thly with the coll i sion velocity, ac cording to the- 
oretic al (e.g.jHarel et al.lll998t iLee et al. 20041) and la boratory 
(e.g., iBeiersdorfer et al.1 [200'lFlMawhorter et al.l 120071) results. 
Nevertheless, the precise dependence on the collision veloc- 
ity is still rather uncertain for a number of ions, and in gen- 
eral it is not expected to vary significantly for speeds above 
400 km s _1 , according to the same works. Moreover, prelimi- 
nary tests with velocity-dependent cross sections revealed neg- 
ligi ble differences in o ur simulation results, in accordance with 
the lGunell et al.l d2004l) conclusions on the matter as well. 

3. Simulation results 

We calculated CX-induced X-ray emission maps in the same ge- 
ometry conditions as the XMM-Newton 2003 observation for 
the most abundant S W ions (listed in Table [TJ, in particular the 
ones whose emission lines were identified in the XMM data. On 
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this date, Mars was at a helioecliptic longitude of about 27° and 
the Earth (XMM) at a longitude of 57°, at a phase angle (Sun- 
Mars-Earth) of -40°. 

3.1. Emission morphology 

In the top panel of Figure|4]we present the projected map of the 
total power flux (W irT 2 sr ) for the spectral lines included be- 
tween 0.3 and 0.9 keV. The white curve superimposed on the 
map represents the average position of the MBS in the same 
projection as the XMM observa tions as inferred fr om the Mars 
Global Surveyor measurements dVignes et al.ll2000T) . In the mid- 
dle panel of Figure |4] we compare the size of our simulation box 
to the field-of-view of RGS that was used to extract X-ray spec- 
tra from Mars . In the bottom panel we reproduce figure 8 from 
iDennerl et al.l d2006l) . showing the total (fluorescence in orange 
and SWCX in green and blue) emission from Mars, for compar- 
ison. 

The model emission displays the characteristic crescent-like 
structure on the sunward side of the planet, where the plasma 
discontinuity at the MBS and the flow of SW around the planet- 
obstacle are clearly identified. Indeed, as shown in the top panel 
of Figure|U the discontinuity in the predicted SWCX emission in 
the model is consistent with the average po sition of the MBS as 
measured by the MGS (IVignes et al.ll2000l) . However, while the 
MGS-inferred MBS is axisymmetric around the Sun-Mars axis, 
the crescent-like structure of the modeled emission is asymmet- 
ric, owing to the influence of the motional electric field. The 
MBS boundary in the simulations (as illuminated by the SWCX 
emission) is pushed towards the Sun in the southern hemisphere, 
and pushed towards the planet near its north pole. More details 
on the plas ma boundaries and the outp ut of the hybrid code can 
be found in lModolo etal] (I2005ll2006l) . 

The overall shape of the model emissio n generally looks sim - 
ilar to the observed total emission pattern (IDennerl et al.ll2Q06b . 
The brightest part of the emission is observed by projection on 
the planet's disk (again, Mars is observed at a 40° phase angle), 
and seems to be produced very close to the planet's surface at the 
limb in the sunward direction. The fainter halo is observed up to 
4 Mar s radii in the north-south direction, while IDennerl et alJ 
claim a much larger extend, up to 8 Mars radii. Due to 
the asymmetry of the MBS, the emission pattern is brighter and 
more elongated toward the south pole of the planet, along the di- 
rection opposite the electric field vector, while there seems to be 
a slight 'depression' of the emission towards the northern hemi- 
sphere of the planet. 

All the individual ion emission maps (not presented here) 
show an equivalent emission pattern. Although the strong asym- 
metry predicted by the model could have been perceived as a 
lobe-like structure by RGS considering the detector's spatial res- 
olution (similar to th e oxygen line maps presented in Figure 7 of 
IDennerl et ail I2006I) . overall the emission pattern seems rather 
continuous in the simulations. The qualitative agreement be- 
tween the emission pattern in the simulations, and the data seems 
to support the idea that the observed tilt of the upper and lower 
wings in the data inde ed has a morphologi cal nature (see discus- 
sion in Section 5. 1 of IDenner l et al. 20Qa). 

3.2. Luminosity and spectral analysis 

The list of lines and their photon fluxes calculated with the model 
are presented in Table[2] In column 1 we list the emitting ion, and 
in columns 2 and 3 the identified transition and energy (in eV) 
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Fig. 4. Top - Total CX-induced X-ray power flux map from Mars 
in units of W irT 2 sr , in the same projection as the XMM- 
Newton observation, at a <p — 40° phase angle. Axes are in Mars 
radii. The planet is represented by the white circle, and the white 
curve marks the MBS position based on MGS data JVignes et al.l 
|200Q|) . The B-Vsw plane corresponds to the ecliptic and the elec- 
tric field (E) axis points north. North p oints downward for direct 
comparison with IDennerl et all d2006i)'s figure 8 reproduced in 
the bottom panel. Middle - Size of our simulation box projec- 
tion at a 40° angle (dashed square), compared to Mars' X-ray 
disk (r < 10") and halo (15" < r < 50") sizes defined by 
IDennerl et al.l d2006l) . Mar s is represented by th e dotted circle. 
Bottom - Reproduction of IDennerl et al.l d2006h 's figure 8. The 
scale is not the same as in the top two panels. (A color version 
of this figure is available in the online journal.) 
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respectively. Columns 4, 5, and 6 list the simulated line fluxes for 
the complet e simulation box, disk , and halo. We adopt the same 
definition as iDennerl et al.l d2006l) for the disk (r < 10") and the 
lower altitude limit for the halo (15" < r), while the upper limit 
for our halo is defined by the simulation box limits (see Figure[4}- 
bottom). Here we compare only CX emission lines; i.e., fluores- 
cence lines are not included in the analysis. Model photon fluxes 
are computed separately for the disk and halo components, as 
well as for the total box size and given in units of 1CT 6 cirT 2 s _1 . 
Finally columns 7, 8, and 9 list the fitted photon fluxes in the CX 
identified lines from the XMM-Ne wton observation (extracted 
from Table 2 in IDennerl et al]|2006l) for the total, disk, and halo 
regions. 

Obviously, although the model includes all the major tran- 
sitions detected by XMM, it widely underestimates the data 
line fluxes, which are in general more than one order of mag- 
nitude higher than the model line fluxes. For energies between 
0.3 and 0.9 keV, the total luminosity within the simulation box 
is 0.23 MW, while for the disk we find 0.06 MW and for the halo 
0.12 MW. 



3.3. The OVII triplet 

The ratio G = (OVII-f +OVII-i)/OVII-r of triplet-to-singlet tran- 
sitions of the He-like OVII ion is used to identify the emission 
mechanism and is usually less than one for hot plasmas where 
excitation of +6 ions occurs due to electron i mpact. This ratio is 
in gen eral greater than three in the case of CX dKharchenko et al.1 
120031) . In the XMM 2003 observations of Mars, the ratio is found 
G~ 6 for the total emission and G ~ 5 for the halo region. Our 
model reflects the atomic data we have used, and yields a mean 
ratio of ~ 8 in t he simulation box, w hich is consistent with the 
values found by_ Dennerl et al.l (12006b . 

However. IDennerl et al.l (120061) report a much fainter G ratio 
(0.8 ±0.6) for oxygen lines from the disk. This interesting result 
may be interpreted in two ways: either a different mechanism 
is responsible for the disk emission or the G factor differs from 
its theoretical nominal value for some physical reason. Because 
the forbidden line emission disappears, it suggests an effect of 
collisional depopulation. We did investigate this possibility. 

If N,„ is the density of the ions in the metastable level, then 
the combined effects of production through CX collisions and 
depopulation by radiative decay and collisions imply that 



For the steady state solution, dN m / dt = 0, 



dN m 1 + U) qC Ty 

- <T cx nVNj - N m 



dt 



(7) 



where cr cx is the CX cross-section, n the density of neutrals (H, 
O), V is the collision velocity, N, the ion density (0 +v in this 
case), a> qc represents the quenching collision frequency, which 
includes all types of collisions and all target atoms, electrons 
and ions able to induce the quenching, and t, is the lifetime of 
the metastable state against radiative decay. 

Quenching collisions remove excited electrons from the 
long-lived metastable states and, because of that, suppress their 
emission. Degradation of the X-ray emission of metastable ions 
may be caused by subsequent CX collisions of metastable +6 * 
ions with atmospheric gas. These CX collisions lead to a forma- 
tion of doubly excited ion states +5 ** that decay mostly due to 
Auger process. If cr q is a cross-section representative of all colli- 
sions for all colliding particles i of total density n q in the planet's 
atmosphere, then, 



Jqc 



Z» 



N m 



a cx nVNi 



1 + OJ qc T r 

and the line radiation flux per unit volume is 
h 



1 + U) qc T r 



(9) 



(10) 



where I is the intensity in absence of collisions. 

Collisions become important and significantly reduce the 



radiation when u) qc j r 



1, i.e. n fl 



(o- q VT r y l . For the life- 



qi" q, V = n q cr q V . 



(8) 



time of the +6 * metastable state t> = 10 3 s, ion velocities 
through the atmosphere of V = 400 km s _1 , and assuming cr q 
= 10~ 15 cm 2 based on a simple geometrical scaling and cross- 
section values of the CX collisions of +6 ions with the CO2 
and CO (5. 10~ 15 cm 2 ), the corresponding critical density is n q 
= 2.5 10 10 cirT 3 . This critical density is relevant to the altitude 
of 150 km where the atmospheric gas is mostly represented by 
CO2 molecules. The same calculation for ions fully decelerated 
down to 4 km s _1 gives n q = 2.5 10 12 cirT 3 , which corresponds to 
altitudes on the order of 100km. The cross-section cr q slightly in- 
creases at low ion velocities, but it would not change the quench- 
ing altitude estimations dramatically because of the exponential 
character of the atmospheric density profiles. 

As a conclusion, quenching of the 2 3 Si OVII metastable 
state and degradation of the forbidden line OVII-f, i.e. relative 
increase in the intensities of the other lines, such as the OVII-r 
resonance line, may happen at altitudes on the order of or be- 
low 150 km, below the Martian exobase. Those regions should 
be able to generate oxygen triplet G values much lower than the 
low density value G~6. (Similar calculations for the quenching 
of the +6 metastable state by ionospheric electrons show that 
this phenomenon is negligible). 

Interpretation of the XMM-Newton results in terms of 
quenching effect s requires that the m ain part of the disk emis- 
sion as defined in lDennerl et al.ld2006l) is generated at an altitude 
of, or below, 150 km. Although our simulations do not include 
the quenching effect, we can estimate where the bulk of the disk 
emission comes from. In the <p=40° projection presented here, 
this is not a trivial calculation, since some of the brightest spots 
on the planet's disk are actually a sum of the line-of-sight (LOS) 
through several layers of the Martian atmosphere. Therefore, we 
estimated the contribution of the layer 0-150 km at the limb 
of the planet in a projection perpendicular to the SW velocity 
vector, knowing that along this direction, the LOS is tangent to 
the most emissive parts of the Martian atmosphere in the subso- 
lar point. With this approach we found that the 0-150 km layer 
contributes no more than 15% of the total emission between 0- 
2180 km (the 10" radius limit of the disk region). However, the 
critical density altitude of 150 km is close to our simulation grid 
resolution, and therefore specific simulations with finer resolu- 
tion, including the quenching effect, are needed to study this hy- 
pothesis in detail. 



4. Discussion and conclusions 

We have modeled the SWCX X-ray emission produced in the 
Martian exosphere in SW conditions very close to the ones 
occurring during the 2003 XMM-Newton observation of Mars 
(IDennerl et al.ll2006h . We calculated projected maps of the to- 
tal SWCX power flux between 0.3 and 0.9 keV, which clearly 
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Table 2. SWCX emission lines in the simulations and RGS spectra of Mars and its halo 



Ion Transition" Energy (eV) Line Photon flux (10 6 cm 1 s ') 

Model Data 6 









Total 


Disk 


Halo 


Total 


Disk 


Halo 


NeVIII 


2p ->ls c 


872.5 








(1.8±1.8) 


(0.5+0.5) 


3.8+2.2 


OVIII 


7p — >ls 


853 


4.2e-6 


1.0e-6 


2.3e-6 








OVIII 


6p — >ls 


846.6 


5.6e-5 


1.4e-4 


3.1e-4 








OVIII 


5p — »ls 


835.9 


1.9e-3 


4.7e-4 


l.le-3 








OVIII 


4p — »ls 


816.3 


6.3e-4 


1.5e-4 


3.5e-4 


3.9+1.5 


(0.4+0.4) 


(1.6+1.6) 


OVIII 


3p — >ls 


774 


1.7e-3 


4.2e-4 


9.6e-4 








OVII 


6p — >ls 


720.9 


9.0e-5 


2.8e-5 


4.3e-5 








OVII 


5p — >ls 


712.3 


2.3e-5 


7.0e-6 


l.le-5 


7.9+2.1 


(0.4+0.4) 


6.2+1.7 


OVII 


4p — >ls 


697.8 


0.01 


4.3e-3 


6.8e-3 








OVII 


4p (T) -^ls 


697 


1.3e-4 


4.0e-5 


6.3e-5 








OVII 


3p — >ls 


665.6 


0.02 


5.9e-3 


9.3e-3 








OVII 


3p (T)->ls 


663.9 


2.4e-4 


7.5e-5 


1.2e-4 








OVIII 


2p — »ls 


653.1 


0.01 


3.3e-3 


7.5e-3 


9.9±2.4 


4.3+2.1 


7.4+2.1 


NVII 


6p — >ls 


648.2 


4.9e-4 


1.3e-4 


2.6e-4 








NVII 


5p — »ls 


640 


6.0e-3 


1.5e-3 


3.2e-3 








NVII 


4p -»ls 


625 


3.0e-3 


7.7e-4 


1.6e-3 








NVII 


3p — »ls 


592.6 


5.2e-3 


1.3e-3 


2.8e-3 








OVII-r 


2p —>ls 


574 


0.03 


7.5e-3 


0.01 


(2.3±2.3) 


7+3.3 


(1.4+1.4) 


OVII-i 


2p (T)->ls 


568.5 


0.03 


0.01 


0.02 


1.5+1.5 


(1.2+1.2) 


2.2+2.2 


OVII-f 


2p (T)->ls 


560.9 


0.21 


0.06 


0.10 


12.8+3.2 


4.2+2.9 


5.3+2.3 


NVI 


4p ->ls 


521.6 


0.02 


5.6e-3 


0.01 








NVI 


4p(T)^ls 


520.9 


8.6e-5 


2.4e-5 


4.5e-5 








NVII 


2p — >ls 


500 


0.05 


0.01 


0.02 


9.3+2.7 


4.5+2.2 


7.1+2.2 


NVI 


3p — >ls 


497.9 


0.01 


3.9e-3 


7.2e-3 








NVI 


3p (T)-»ls 


496.5 


7.8e-5 


2.2e-5 


4.1e-5 








CVI 


6p — *ls 


476.2 


3.7e-4 


9.5e-5 


2.0e-4 








CVI 


5p — >ls 


470.2 


0.03 


8.4e-3 


0.02 


5.5+2.2 


2.7+2.3 


5.5+1.8 


CVI 


4p -»ls 


459.2 


0.16 


0.04 


0.08 


4.5+3.2 


5.7+2.8 


5.3+2.5 


CVI 


3p — »ls 


435.4 


0.10 


0.03 


0.05 


6.6+2.3 


(1.9+1.9) 


7.9+2.7 


NVI 


2p — >ls 


430.7 


0.04 


0.01 


0.02 








NVI 


2p(T)^ls 


426.1 


0.04 


0.01 


0.02 








NVI 


2s (T) ->ls 


419.8 


0.27 


0.08 


0.14 








CV 


5p — »ls 


378.5 








9.9+4.2 


(2.8+2.8) 


8.4+3.1 


cv 


4p ->ls 


370.9 


0.03 


0.01 


0.02 








CV 


4p(T)^ls 


370.4 


5.5e-5 


1.8e-5 


2.8e-5 








CVI 


2p — >ls 


367.4 


0.61 


0.16 


0.33 


16.2+5.2 


11+4.2 


9.1+3.5 


cv 


3p — >ls 


354.5 


0.08 


0.02 


0.04 








cv 


3p (T)-»ls 


353.4 


1.8e-4 


5.7e-5 


9.0e-5 








OVIII 


2s — >ls 


326.5 


1.4e-3 


3.3e-4 


7.6e-4 








cv 


2p — >ls 


307.9 


0.16 


0.05 


0.08 








cv 


2p (T)->ls 


304.2 


0.07 


0.02 


0.04 








Total luminosity (MW) 




0.23 


0.06 


0.12 


11. 8 C 


5.8'' 


9A d 



Notes. {a) Transitions involving a triplet state are n oted with a (T) . m On ly the SWCX emission lines from lDennerl et af] d2006l) are reported in 
this table. w This line has been falsely identified bv lDennerl et al] d2006h as the NeVIII 2p — >ls transition because it would involve the electron 
being captured in an excited core level, which cannot be produced by CX. (< " Excluding the NVII 3p— >ls due to its large uncertainty. 



show the bow shock formed at the SW plasma - planetary neu- 
tral interaction region. Just as in mo deling of the SWCX e mis- 
sion from the Earth's magnetosheath dRobertson et al.l2006l) . our 
model maps demonstrate the potential of SWCX X-ray emission 
in the studies of electromagnetic structures and SW-neutral in- 
teractions around solar syste m bodies, in support of proposed 
missions of X-ra y imagers dBranduardi-Ravmont et al.l 1201 it 
iKuntz et al.ll2Q08b . 

We find a total luminosity of 0.23 MW in the energy range 
between 0.3 and 0.9 keV, almost two orders of magnitude 
lower than the XMM-Newton data (11.8 MW), but closer to the 
Chandra ob servations of th e planet (0.5+0.2 MW in the range 
0.5-1.2 ke V. |Pennerlll2002l) . The bulk of the emission in our re- 
sults is found at distances up to 4 Rm, compared to 3 and 8 Rm 
for the Chandra and XMM-Newton data, respectively. 



Previous simulations from iHolmstrom et al.l (l200lh with 
similar SW conditions (V m = 400 km s _1 , n m = 2.5 cm -3 ) 
yielded 2.4 MW (1.52 10 25 eV/s) for solar minimum and 
1.5 MW (0.93 10 25 eV/s) for solar maximum conditions of the 
neutral coro n as, w ithin a 10 Rm radius. On the other hand, 
iGunell et all d2004l) find 1.8 MW for a 6 R M simulation box 
(with V sw = 330 km s _1 , n JVV = 4.4 ctrT 3 inputs). Our simula- 
tion box only extends up to a 43" radius (~7 Rm, see bottom of 
Figure|4]i, but a crude upper limit estimate of the additional emis- 
sion in the annulus between 43" and 50" (~8 Rm) only yields an 
additional 5% compared to the total luminosity in the box. 

Besides the small differences in the proton fluxes used and 
the size of the simulation boxes, two main effects can account 
for the divergence of our results from previous simulations. First, 
all previous simulations have assumed average slow SW abun- 
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dances for the minor ions. These abundances are from three to 
five times higher than the real-time data we used in our simula- 
tions (see Table Q]) and should have a significant effect (approxi- 
mately of the same factor) in the comparison with previous simu- 
lations. If we apply an average abundance three times higher for 
all the ions included in our simulations, then the total luminosity 
reaches ~0.7 MW in the simulation box. 

A second effect that may, in part, explain the discrepan- 
cies with previous simulations is the lack of the H2 a nd He 
cor onas in our simulations. A ccording to iFox & Had (120091) 
and iRrasnopolskvl d2002l 1201 Oh . the H2 and He density in the 
Martian atmosphere can be comparable to the H density, de- 
pending on solar conditions. For solar maximum, the H2 density 
is almost an order of magnitude higher than H at 400 km alti- 
tude (but still an order of magnitude lower than O), while He 
can be anywhere between two and ten time s the H density at 
the same altitude (s ee profiles extracted from[Fox & Hac 20091: 
lKrasnopolsky||201d in Figure[3]). Since the H2 and He corona is 
not included in our simulation, this may have an important effect 
on the SWCX luminos ity, especially in the disk region. From 
iGreenwood et al.l d2001l) we deduce that the CX cross-section of 
highly charged C, N, and O ions with H2 is near to the cross- 
section for H. Therefore, based on the crude assumption that 
the SWCX emission is proportional to the neutral density, and 
had we included the H2 and He coronas in our simulations, the 
SWCX luminosity could be, at best, two or three times as strong 
in the disk (assuming that the H2 and He sum to 20 times the H 
density in the disk). However, the H2 and He coronas are domi- 
nant only at low altitudes, up to about 1000-2000 km. After that, 
the H2 and He densities decrease more rapidly than the H den- 
sity, and have much less effect on the halo emission. In future 
work, we intend to include the H2 and He coronas and perform 
a parametric study of the solar activity influence on the SWCX 
emission from Mars. 

If we sum up the estimated additional emission from the 
missing H2 and He coronas in the disk and the 5% of the outer 
halo region between 7 and 8 Rm, the total luminosity in our sim- 
ulations could be on the order of 0.4 MW under the same SW 
conditions. If we also adjust to the average slow SW ion abun- 
dances, then our simulated luminosity can reach between 1 .2 and 
2.0 MW, which bri ngs us to almost perfect agreement with pre- 
vious simulations (lKrasnopolskvll2000l: iHoimstrom et alJl200lt 
iGunell et alJl2004h . but still a far way from the observed values. 

Another interesting hypothesis for the extended SWCX halo 
around Mars observed by XMM could be the presence of a hot 
population of neutral atoms at large dist ances from the planet. 
Indeed, according to iDelva et all d201 ll) . upstream proton cy- 
clotron waves have been observed by up to 11 Rm- These waves 
are generated by pick-up of planetary neutrals by the solar wind, 
and therefore their occurrence at large distances may come from 
the presence of an extended hydrogen exosphere around the 
planet. The large extent of this hot neutral component is fur- 
ther supported by the fact that the excitation of these waves to a 
measurable amplitude requires an accumulation of ions (gener- 
ated by ionization of planetary neutrals) in a plasma cell, while 
it is convected by the S W. Therefore these ions need to be gener- 
ated by a large quantity of neutrals located upstream in the SW 
with respect to the position where the cyclotron waves are ac- 
tually observed. The influence of this hot neutral corona to the 
ion feedback and electromagnetic interface around the planet has 
not yet been included in any model of the SW - neutral interac- 
tion of the Martian environment to our knowledge, and should 
be explored in the future. 



CORHEL simulations at Mars 




16/11/03 18/11/03 20/11/03 22/11/03 24/11/03 
Date 



Fig. 5. Top - Simulated proton density (solid) and velocity 
(dashed) calculated with the CORHEL model for the Mars or- 
bit for the2003/l 1/18 CME propagating through interplanetary 
space. According to the simulations the perturbed CME plasma 
arrived at Mars about a day and a half after the XMM-Newton 
exposure (represented by the vertical lines). Bottom - Same as 
top only calculated for the Earth orbit. ACE 1-day gliding av- 
erages of density and velocity measurements are also plotted in 
gray for comparison. 



An additional explanation for the discrepancies between the 
simulation results and the XMM-Newton data is the large uncer- 
tainty of extrapolating Ll -based in-situ SW measurements to the 
Mars orbit. If we consider the longitudinal difference A6 ~ 30° 
in the positions of Earth and Mars, then the time delay due to 
the radial propagation of S W particles is somewhat compensated 
for by the ~27-day solar rotation (At = y- — 13 3^ flv X since 
Mars is behind the Earth during this period. In that case, the 
ACE SW window is situated approximately between 2003/1 1/20 
14:00 UT and 2003/11/21 06:00 UT. In this interval the ACE 
data show a strong density and abundance enhancement related 
to a halo coronal mass ejection (CME) emitted on 2003/11/18. 
The SW ion flux N, measured in this interval was on average 
18 times higher than the real-time values we used in our simula- 
tions, which could yield a total luminosity of 4 MW (~6.3 MW, 
if we include the H2 and He coronas as well), in better agreement 
with, but still quite lower than, the observed ~12 MW. 

However, a different proof argues against this hypothesis. In 
Figure [5] we present results from a CORHEL simulation rele- 
vant to the 2003/1 1/18 CME propagation in interplanetary space. 
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The simulation results were extracted from the Community 
Coordinated Modeling Center (CCMCQ) Runs-on-Request sys- 
tem. The top panel of Figure [5] presents the simulation results 
for the CME impact on Mars, and we clearly see that the CME 
arrived a day and a half later than the XMM-Newton exposure 
window. 

To test the timing accuracy of the CORHEL simulations we 
also checked the simulation results on Earth's orbit and com- 
pared them to the ACE data (Figure[5]-bottom). Although the ab- 
solute velocity values in the simulations have up to 200 km s _1 
difference compared to the actual data measured by ACE at the 
LI point, the timing seems accurate, and so it seems to indicate 
that the extrapolation to Mars is also accurate and that the CME 
missed the XMM-Newton observation window at Mars indeed. 
Even if we apply the data-model velocity difference to estimate 
the error-margin of the CORHEL simulation, the CME reaches 
Mars just at the end of the XMM exposure, still too late to ef- 
fectively influence the SWCX emissivity of the planet. Based 
on these tests, and owing to the lack of in-situ plasma measure- 
ments at Mars simultaneous to the XMM-Newton observations, 
it is impossible to assert the actual impact of the CME on the 
planet's SWCX emission. 

For the various reasons discussed above, we feel that the 
controversy over the very strong SWCX emission detected at 
Mars with XMM-Newton still remains, especially since it has 
not been verified by subsequent X-ray observations of the planet. 
However, all subsequent observations of the planet until now 
were scheduled during a peculiar solar minimum, when solar ac- 
tivity and SW flux were particularly low. The issue should hope- 
fully be elucidated if a similar event is reproduced. The condi- 
tions for such an event are becoming favorable at present and as 
we approach the solar maximum in the next couple of years, so 
new X-ray observations of Mars are needed in order to constrain 
the SWCX emission of the planet. 

The discrepancy between the model and the XMM data can- 
not be lifted through the present analysis, but the sensitivity of 
the model results to the neutral corona parameters shows that X- 
ray observations combined with simulations allow the Martian 
and planetary exospheres, in general, to be studied remotely. 
This can potentially have interesting applications for studying 
extrasolar planet exospheres and their mass-loss processes. 

Several recent studies report the presence of extended hy- 
drogen atmospheres (as well as C and O, in some cases) 
around extrasolar planets, revealed by means of absorption of 
the stellar H I Lyman-a (C II and O I, respectively, for the 
other compo nents) line when the planet transits in front of 
the star (e.g., IVidal-Madiar et al.|l200l lEhrenreich et al.ll200l 
iLecavelier des Etangs et al.l 2010 ). An interesting debate arose 
following the H detection, with the first interpretation sup- 
porting the idea that H atoms are escaping the planet's ex- 
osphere, owing to hydrodyn amic blow-off, and acce l erated 
by stellar radiation pressure (fVidal-Madiar et all 120031 l2004t 
ILecavelier des Etangs et al.ll2008l) . A different approach i ntro- 
duced bv lHolmstrom et al.l (120081) . and further explored bv lGuol 
d201 lb and lLammer et al.l(l201 ll) . proposed that radiation pres- 
sure alone cannot suffice and the Lyman-a absorption could be 
explained by energetic neutral atoms (ENAs) produced in charge 
exchange collisions between the exospheric neutrals and the stel- 
lar wind protons. In this case, ENAs would allow the study of 
the plasma properties in the stellar wind, and while they do 
not exclude a large atmospheric escape, they cannot provide 



http://ccmc.gsfc.nasa.gov 



any strong constraint on it (see reply by Holmstrom et al. to 
ILecavelier des Etangs et al"1l2008h . 

If the star is endowed with a solar-like wind, then CX colli- 
sions are bound to occur between the ionized wind particles and 
exospheric neutrals. In that case, CX-induced X-rays will also 
be produced if the wind's high ion content is at least solar-like. 
Therefore, in principle, it would be possible to study the interac- 
tion of the stellar wind with the planet's exosphere by CX X-ray 
imaging and spectroscopy. This would provide a complementary 
method, less ambiguous than Lyman-o' absorption, to study both 
the plasma and neutral properties of the planet's environment. 
In theory, the CX X-ray spectroscopy would help to put con- 
straints to the stellar wind ion mass flux and neutral composition 
and density of the ex osphere, since the CX lin e ratios depend on 
the collision energy dKharchenko et al.ll2003l) and the CX cross 
sections depend on the neutral target dGreenwood et al.ll2.00lb . 
In addition, CX X-ray imaging could constrain the magnetic 
boundaries around the planet, giving information on whether the 
planet is strongly magnetized or not. 

Although the study of star-planet interaction with CX- 
induced X-ray imaging and spectroscopy seems promising, we 
feel that technically it is still unfeasible, as it requires instru- 
ments with large collecting areas, and enough spatial and spec- 
tral re solution to separate the planet's X-ray signal from the 
star's. IWargelin & Drake! (1200 ll) advanced this CX X-ray based 
approach to study the mass-loss rate of late-type dwarfs through 
the CX interaction of their winds with neutral material of the IS 
medium (in the same context as heliospheric SWCX emission), 
but did not detect any significant CX halo in t he cases of Proxima 
Centauri and the M dwarf R oss 154 studied (IWargelin & Drake] 
12002b IWargelin et al.1120081) . Recent studies of star-planet inter- 
actions in X-rays have mainly focused on the star's X-ray envi- 
ronment to determine the irradiation-driven mass loss of super- 
Earths or the effe cts of giant planets on the chrom ospheric and 
coronal activity dPoppenhaeger et all 1201 ll 1201 2|) . concluding 
that "observational evidence of star-planet interaction in X-rays 
remains challenging". Hopefully, the next-generation X-ray ob- 
servatories proposed (e.g., ATHENA) will be able to shed more 
light in star-planet interactions in extrasolar systems. 
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